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We present an indirect imaging method that measures both amphtude and phase information 
from a transmissive target. Our method is based on an optical eigenmode decomposition of the 
light intensity and the first-order cross correlation between a target field and these eigenmodes. 
We demonstrate that such optical eigenmode imaging does not need any a priori knowledge of the 
imaging system and corresponds to a compressive full-field sampling leading to high image extraction 
efficiencies. Finally, we discuss the implications with respect to second-order correlation imaging. 
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The last two decades have seen the emergence of in- 
terest in schemes for interaction-free imaging [1] and in- 
direct imaging [2] in both the quantum and classical do- 
main [3]. Such schemes inherently use light fields to im- 
age that have themselves never scattered from the object. 
Prominent amongst these schemes has been the concept 
of ghost imaging [4-7]. These methods allow the image 
of an unknown object to be nonlocally reconstructed by 
the intensity correlation measurements between two light 
fields. The light that illuminates the object is typically 
collected by a single pixel detector that itself has no spa- 
tial resolution. The underlying physics of ghost imaging 
has seen major debate about its classical and quantum 
implications as both entangled source and thermal light 
can be used. Whilst a powerful concept, ghost imaging 
does not inherently reveal phase information about an 
object as it relies on a second-order correlation for its 
implementation [2, 8]. Additionally, the resolution of the 
recovered images is determined by the size of the speckle 
placed on the object plane [4]. Any imperfections or aber- 
rations within the system may not be readily dealt with 
in such a system. A significant step forward would be an 
indirect imaging scheme based on correlation measure- 
ments that could retrieve both amplitude and phase and 
inherently accommodate issues relating to aberrations or 
imperfections within the optical system 

In this paper, we realise such a scheme using the con- 
cept of optical eigenmodes [9, 10]. We split the laser light 
in two different beams. One beam does not interact with 
the target, but illuminates a high-resolution CCD cam- 
era (multi-pixel detector). The other one interrogates, in 
transmission, the target (or sample) and then illuminates 
a photodiode (single-pixel detector) providing no spatial 
resolution. The transmission wave-front of this beam is 
decomposed, using an optical lock-in amplification tech- 
nique, onto an orthogonal set of optical eigenmodes. The 
lock-in amplification corresponds to performing a first or- 
der cross-correlation and as such is distinct from present 
ghost-imaging techniques. In turn this leads to the re- 
tention of the phase information of the object. Further, 
our approach foregoes the point-by-point scanning and 
allows a rapid full field image extraction [11]. 



The outline of the paper is as follows. Firstly, the con- 
cept of the optical eigenmodes is briefly reviewed and 
then we discuss the problem of combining eigenmodes 
and indirect imaging. We introduce the experiment and 
show an application of the eigenmode method for indi- 
rect imaging of a target. We demonstrate the advantage 
of optical eigenmode imaging in terms of resolution and 
phase information. Finally, we show the relationship be- 
tween optical eigenmodes and the second-order correla- 
tion function. 

Method. To deflne the intensity optical eigen- 
modes [10], we decompose a linearly polarised electro- 
magnetic fleld E into a superposition of N monochro- 
matic (e^^*) "test" fields: 

E = 5^a*E,- ; E* = ^E^a,. (1) 

j k 

Here, we consider the field intensity m^^^ integrated over 
a region of interest (ROI), as defined by: 

m^i)^E)= [ (icT E-E* (2) 

valid for linearly polarised light and not tightly- focused 
beams. The ROI represents the detector active area. 
Equation (2) can be written in a general quadratic matrix 
form: 

m(^)(E)=^a*M,fcafc (3) 

where the elements Mjj. are constructed by combining 
the fields Ej and Ek for j, k = l...A^: 

Mjk= [ daEj.El. (4) 

-'ROI 

The optical eigenmodes are defined by: 

^^ = ^J2^h^^ ; ^^* = ;;^E^«E* (5) 

with 

^ MjkVij = X^Vik (6) 
j 
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FIG. 1. The experimental set-up used to perform indirect 
imaging with the optical eigenmode formalism. The inset 
shows a scheme of the target of three holes used in our exper- 
iment. Abbreviations: L=lens, M=mirror, BS=beam splitter, 
FP= focal plane, PD= photodiode. 



where is an eigenvalue and V£j the associated eigen- 
vector. The matrix Mjk shows two important proper- 
ties. Firstly, it is Hermitian meaning that all the eigen- 
values (A^) are real and can be ordered where A^^^ is 
the largest eigenvalue, A^^^ the second largest eigen- 
value, etc... . Secondly, two eigenvectors correspond- 
ing to different eigenvalues are orthogonal which means 
that: /j^oi ' ^fc ~ ^jf^' matter of fact, each of 

these eigenvectors V£j corresponds to a superposition of 
the initially considered fields. An unknown field T can 
be decomposed onto the eigenmodes using its projection 
defined by: 
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daT- 



(7) 



where q corresponds to the complex decomposition co- 
efficients of the field T in base E^. If the fields form a 
complete base, we can perfectly reconstruct the unknown 
field T from the projection using T = c|E^. We remark 
that the completeness of the base is dependent on the 
initial fields probing all the degrees of freedom available. 

Optical eigenmode imaging is based upon a direct mea- 
sure of the projection coefficients c| and the experimental 
reconstruction of the unknown field T, that is essentially 
the transmission function through the target. 

Experimental. The basic set-up used for the optical 
eigenmode imaging is shown in Fig. 1. A He-Ne laser 
(Thorlabs, A = 633 nm, Fmax = 10 mW) beam is ex- 
panded to fill the aperture of the spatial light modulator, 
SLM (Hamamatsu LCOS-SLM X10468, resolution = 600 
X 800 pixels, pixel size = 20 x 20 /im^. Refresh rate = 60 
Hz) and then phase-modulated in the first-order configu- 
ration [12]. This configuration allow us to achieve phase 
and intensity modulation using the interfering efficiency. 
The phase-modulated beam is divided in two symmetric 
parts by a 50:50 beam splitter. The reflected beam inter- 
rogates a transmissive target placed in the focal plane of 



the Fourier lens (L3=75cm). A second Fourier lens (L4) 
integrates the transmitted light on a photodiode provid- 
ing no spatial resolution. The transmitted beam is fo- 
cused directly on a high-resolution CCD camera (Easier 
pilot piA640-210mg, resolution = 648 x 488 pixels, pixel 
size = 7.4 X 7.4 /im^) without interacting with the tar- 
get. By projecting the optical eigenmodes on the target 
it is possible to reproduce the target image on the CCD 
camera. 

The experiment consists of two different steps. Firstly, 
it is necessary to determine the optical eigenmodes (E^) 
in the CCD camera plane by interfering the test fields 
Ej created in the SLM plane. Next, we select the n 
most intense optical eigenmodes and interfere them with 
a reference signal (Eref). The reference signal (Eref) is 
random phase encoded [13] on the same SLM and thus 
shares the same optical path as the test field. In this 
way, we can reconstruct the amplitude and phase of the 
linearly polarised fields Fj on the CCD camera: 



R-l 



e-'2-p/%.|2 



(8) 



p=0 



where R is an integer number greater then 3. This pro- 
cedure can be seen as a phase sensitive lock-in technique 
where the reference beam E^ef corresponds to a reference 
signal with respect to which the phase and amplitude of 
Ej is measured. This is done by measuring R times the 
intensity on the CCD camera while changing the relative 
phase between the reference and the test field. This in- 
tensities are then combined in a equation (8) in an inverse 
Fourier transform delivering a measure proportional to 
the complex test field Ej. This approach generalises the 
four points {R = 4) method presented in [10]. We re- 
mark here that, due to the lock-in technique, using a 
larger number of measures R experimentally delivers a 
better signal to noise ratio. The intensity operator Mjj^ 
is given by: 



jk 



daFi ■ Ft 



(9) 
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where ROI is a region of interest on the CCD camera. 

In the second experimental step, we project the opti- 
cal eigenmode on the target (T) using a first-order cross- 
correlation with the reference beam. The projection co- 
efficients Ck are related to the signal Sk measured by the 
single-pixel detector (PD) as follows: 



Sk 



i27vp/R 



-T 

= / d(jT- E*ef • / dcr T* • E/e 
JROI ^ROI 



d<7T-(E,ef + e-'2-p/fiEfc) 



da T • E*,f 



(10) 
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FIG. 2. Examples of the principal experimental intensity op- 
tical eigenmodes. The resolution of each image is 120 x 120 
pixels and the scale is not shown for simplicity. 

where /j^qi ^ ' -^ref gi^^s the coupling coefficient to 
the target and is essentially a constant value. There- 
fore, from the signal acquired on the PD, we can measure 
the coefficients which makes it possible to reconstruct 
the complex target field. Finally, using a Dirac distribu- 
tion target field, we can deduce the point spread func- 
tion EpsF which determines the resolution of the optical 
eigenmode imaging method: 

EpsF = y2^j [ ^{r - roWjda (11) 

Results and discussion. In our experiment, we used 
N = 1353 masks consisting of 33 horizontal and 41 verti- 
cal independent deflection angles to determine the inten- 
sity optical eigenmodes (E^) in the CCD camera plane. 
After the Fourier lens (L3), these deflections correspond 
to a regular grid scan over an area of 1.2mm by 1.3mm. 
The reference field (Eref) was chosen to be the median 
beam out of the 1353 beams used. The ffist 20 experi- 
mental intensity optical eigenmodes are shown in Fig. 2 
displaying visually an orthogonal behaviour similar to 
TEM modes. 

From the experimental intensity optical eigenmodes, 
we select the ones with an efficiency above 0.1% with 
respect to the maximum possible, giving us typically 
around 140 eigenmodes. This corresponds to a tenfold 
image data reduction when compared to the initial num- 
ber of test fields {N = 1353). Further, we noticed no 
clear increase in the number of suitable eigenmodes when 
increasing N. The quality of this optical compression 
technique can be assessed experimentally by determin- 
ing the PSF (11) showing the loss of image information 
when reconstructing a Dirac distribution field. Fig. 3-a 
shows the experimental PSF obtained by combining 147 
intensity optical eigenmodes representing the resolution 
limit of the system. Further, we numerically project the 
experimental eigenmode fields onto a simulated Laguerre- 
Gaussian target with a unit vortex charge. The theoreti- 




20 40 60 80 100 120 20 40 60 80 100 120 

Pixels Pixels ^ ^ '5 

. 0.5 3 

Theoretical Phase Experimental Intensity ^ 




FIG. 3. Experimental optical eigenmode Imaging of a 
Laguerre- Gauss (LG) beam, (a) Point Spread Function (PSF) 
obtained by combining 147 Intensity OEi. (b) Theoretical in- 
tensity and (c) phase of an decomposed LG beam (Intensity 
scale bar is from (blue) to 27v (red)), (d) Experimental 
reconstruction of the LG beam using the 147 optical eigen- 
modes. 



cal reconstructed field is displayed in Fig. 3-b and Fig. 3-c 
as field intensity and phase respectively. Fig. 3-d shows 
the intensity on the CCD camera when the SLM dis- 
plays the optical eigenmode superposition corresponding 
to T = X]^^!^^- These figures show that optical eigen- 
mode imaging provides phase and intensity information 
as predicted in equation (7). 

For the indirect imaging we used a target consisting of 
three holes each with a diameter of ~ 150 pm (see Fig.l for 
more details). Fig. 4-a shows the conventional transmis- 
sion image of the target obtained by scanning the laser 
spot, using the initial 1353 deflection masks, through the 
holes. In this figure, we plot the signal intensity detected 
by the single-pixel detector (PD) for each angle of the 
ligth scanning. By contrast. Fig. 4-b and Fig. 4-c show 
the indirect image obtained, respectively through numer- 
ical and experimental optical eigenmode superpositions, 
on the CCD camera. The latter was produced using the 
SLM mask reported in Fig. 4-d. Direct comparison of 
the theoretical superposition (Fig. 4-b) and the experi- 
mental one (Fig. 4-c) clearly reveals a good agreement. 
These results suggest that the numerical and the exper- 
imental distributions verify: oc c/^, demonstrating the 
linearity of our optical system. Further, by comparing 
Fig. 4-a and c, we can clearly observe an enhancement 
in the resolution of the optical eigenmode image with re- 
spect to the conventional transmission image using the 
initial scanning "test" fields. 

Remarkably, the optical eigenmodes described above 
simplify the ffist-order correlation functions. Indeed, we 
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FIG. 4. Indirect Imaging of a three holes target, (a) Conven- 
tional transmission image of the target reconstructed from the 
intensity signal collected by the PD as a function of the beam 
displacement in the target plane, (b) Corresponding numer- 
ical indirect intensity Image of the target; (c) Experimental 
indirect optical eigenmode Image; (d) Final mask encoded on 
the SLM. 

have 

G«(^)= / {E{t)E*{t + T))da 
Jroi 

= e"'"^ /Y.a*M,kak\ =T.G%\r) (12) 

\ j,k I jk 

with 

Gf^{r) = / {E,mi{t + r)) da = e— 5,, (13) 
Jroi 

where the symbol < ... > indicates the ensemble averag- 
ing. This relationship shows that the optical eigenmodes 
are independent with respect to the first-order correla- 
tion function i.e. amplitude and phase fluctuations of 
two different eigenmodes are not correlated. This im- 
plies that any random thermal fluctuations of the illumi- 
nation source can be decomposed into independent fluc- 
tuations each corresponding to an eigenmode. We also 
remark that the phase information is maintained in the 
first order correlation process. It is these two properties 
that make the optical eigenmode indirect imaging pos- 
sible. Additionally, illuminating the target directly with 
the eigenmodes does not rely on random thermal fluc- 
tuations to explore their Hilbert space in an inefficient 
way. 

In conclusion, we have shown that optical eigenmode 
imaging allows the indirect complex reconstruction of 
a transmissive target. This method is based on the 
decomposition of the field into intensity optical eigen- 



modes. The decomposition coefficients are directly mea- 
sured by the first-order cross-correlation between the un- 
known target field T and the optical eigenmodes. Su- 
perimposing the optical eigenmodes with the decompo- 
sition coefficients yields the reconstructed target field T. 
Importantly, the complex decomposition coefficients pro- 
vide phase-information potentially allowing the recon- 
struction of optical path length information of the target 
akin to transmission optical tomography. Further, we 
show that the optical eigenmode imaging corresponds to 
a compressive full-field sampling improving the imaging 
speed without loss of details enabling fast spectroscopic 
imaging applications such as Raman imaging. Finally, 
we used the concept of optical eigenmodes to describe 
the first-order and second-order correlation functions. In 
future work we aim to explore the implications of optical 
eigenmodes in the field of quantum optics. 
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